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We deduce the normal state angle-resolved single-particle self-energy T.(9,uj) and the Eliashberg 
function (i.e., the product of the fluctuation spectrum and its coupling to fermions) a 2 F(8,ui) 
for the high temperature superconductor Bi2Sr2CaCu2 0g+5 from the ultra high resolution laser 
angle-resolved photoemission spectroscopy (ARPES). The self-energy £(0,w) at energy uj along 
several cuts normal to the Fermi surface at the tilt angles 6 with respect to the nodal direction in a 
slightly underdoped Bi2Sr2CaCu2 08+5 were extracted by fitting the ARPES momentum distribution 
curves. Then, using the extracted self-energy as the experimental input, the a 2 F(0, uj) is deduced by 
inverting the Eliashberg equation employing the adaptive maximum entropy method. Our principal 
new result is that the Eliashberg function a 2 F(8, lj) collapse for all 6 onto a single function of lj up 
to the upper cut-off energy despite the 6 dependence of the self-energy. The in-plane momentum 
anisotropy is therefore predominantly due to the anisotropic band dispersion effects. The obtained 
Eliashberg function has a small peak at uj ~ 0.05 eV and flattens out above 0.1 eV up to the angle- 
dependent cut-off. It takes the intrinsic cut-off of about 0.4 eV or the energy of the bottom of the 
band with respect to the Fermi energy in the direction 9, whichever is lower. The angle independence 
of the o?F{6 , uj) is consistent only with the fluctuation spectra which have the short correlation 
length on the scale the lattice constant. This implies among others that the antiferromagnetic 
fluctuations may not be underlying physics of the deduced fluctuation spectrum. 



PACS numbers: PACS: 

I. INTRODUCTION 

It is generally agreed that understanding the normal 
state properties is prerequisite to the high temperature 
superconductivity because the pairing instability is a nor- 
mal state Fermi surface instability. A number of stud- 
ies have been conducted to elucidate the normal state 
charge and spin dynamics and momentum anisotropy of 
the high temperature superconductors*^— The angle re- 
solved photoemission spectroscopy (ARPES), owing to 
its unique momentum and energy resolution, has been 
quite powerful in uncovering the in-plane momentum 
anisotropy of the quasi-particle (qp) dynamics of the 
cupratestS Early ARPES measurements^ on optimally 
doped (OP) Bi 2 Sr 2 CaCu20 8+( 5 (Bi2212) and subsequent 
measurements^— along the nodal cut of (0, 0) — (it, ir) 
direction showed the marginal Fermi liquid (MFL)^i be- 
havior of —ImY>(uj) oc uj and found a kink around 0.06 eV 
in the qp dispersion. Kaminski et al. extended these mea- 
surements to off-nodal cuts to investigate the in-plane 
anisotropy of the qp scattering rate -/mS(S,u) around 
the Fermi surfaced They reported that the functional 



form of the scattering rate for under (UD) and optimally 
doped sample can be written as a + buu, where the elastic 
term a(6) is anisotropic in correlation with the pseudo- 
gap, and the inelastic term, 6, is isotropic around the 
Fermi surface; the scattering rates become isotropic for 
the heavily overdoped (OD) samples. On the other hand, 
Chang et al. reported that both the elastic and inelastic 
terms are anisotropic, exhibiting a minimum along the 
nodal direction^ 

The in-plane anisotropy may also be probed by the 
angle dependent magneto-resistance (ADMR)^ or the 
Raman scattering experiments.— For instance, the po- 
lar ADMR measurements in OD Tl 2 Ba 2 Cu0 6+(5 (T12201) 
were analyzed in terms of the transport scattering rate 
Ttr which consists of the isotropic and anisotropic terms 
as 

r tr (0, T) = T {9) + r a sin 2 (2#)T + r 2 T 2 , (1) 

where To(0) is proportional to \/vp{0), i.e., the in-plane 
density of states The anisotropic T linear term, inter- 
estingly, has the same anisotropic form as the uj linear 
contribution to the ImT, that Chang et al. reported, as 
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mentioned above. 

Besides its own significance, the momentum anisotropy 
of the single-particle self-energy is also necessary to 
understand the enormous amount of spectroscopic and 
transport properties of the cuprates. Moreover this self- 
energy is determined by coupling to the fluctuation spec- 
tra which is the essential physical quantity to understand 
in the cuprates. This may be done using the inversion of 
the Eliashberg equations provided the reliability of the 
approximations in the latter is checked consistently. Re- 
cently Schachinger and Carbotte^» performed this inver- 
sion using the maximum entropy method (MEM) 18 from 
the nodal cut ARPES data on Bi2212 of Zhang et alH 
The obtained fluctuation spectra showed a weak peak 
around 0.06 eV and a cut-off near 0.4 eV, which is sim- 
ilar to that obtained from the frequency dependent con- 
ductivity on Bi2212.- This may be anticipated because 
the conductivity, being proportional to the Fermi veloc- 
ity squared, is strongly weighted by contributions from 
near the nodal region. The analysis of the frequency 
dependent conductivity in terms of the self-energy and 
Eliashberg function has been applied in many cases^ir— 

However, the information obtained from the nodal di- 
rection alone is insufficient to characterize the self-energy 
and the momentum dependence of the fluctuation spec- 
tra. In the present work we analyze the momentum 
anisotropy and frequency dependence of the qp self- 
energy £(0,cj) in a slightly under-doped Bi2212 sample 
of T c = 89 K along the cuts of the tilt angles = 0,5, 10, 
15, 20, and 25 degrees with respect to the nodal direction 
from (-7T, 7r) point in the Brillouin zone as shown in Fig. 
[TJ The angle-resolved photoemission measurements have 
been carried out on a VUV laser-based ARPES system^. 
In the experimental setup, there are two angles, 9 and <fi, 
to be controlled. They are then converted to the actual 
tilt angle 9 and the amplitude k± of the inplane wave- 
vector k. The k± in this paper is referred to the distance 
from the (ir, it) point. The thick bars along the cuts of the 
Fig. [1] indicate the range of k± of the collected ARPES 
data as the angle is varied by approximately 30 de- 
grees. The actual tilt angle 9 deviates slightly from the 
control tilde angle 9 as the <f> is varied. The change of 9 is 
small and is disregarded in the analysis to be presented 
here. See the more technical discussion below Eq. ([3]) in 
section II. 

The photon energy of the laser is 6.994 eV with a band- 
width of 0.26 meV. The energy resolution of the electron 
energy analyzer (Scienta R4000) is set at 1 meV, giving 
rise to an overall energy resolution of 1.03 meV which 
is significantly improved from regular 10^15 meV from 
regular synchrotron radiation systems. The angular res- 
olution is ~ 0.3°, corresponding to a momentum resolu- 
tion ^0.004 A -1 at the photon energy of 6.994 eV, more 
than twice improved from 0.009 A -1 at a regular pho- 
ton energy of 21.2 eV for the same angular resolution. 
The slightly-underdoped Bi2212 single crystals with a 
superconducting transition temperature T c —89 K were 
cleaved in situ in vacuum with a base pressure better 




FIG. 1: The Fermi surface of Bi2212 in the first Brillouin 
zone. The blue solid curve is the FS from Eq. and the 
solid dots are the experimentally determined FS at 6 — 0, 5, 
10, 15, 20, and 25 degrees. is the distance from the (tt, 7r) 
point. The thick bars along each cut indicate the ranges of 
experimentally measured ARPES MDC data. Parts of them 
which do not deviate much from the lines were actually used 
in the TB fitting as indicated in Fig. 2] 
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FIG. 2: The temerature dependence of the resistivity of a 
slightly underdoped Bi2212 of T c = 89 K. T* is estimated 
from the deviation from the linearity to be T* ~ 160 K. 



than 5xl0 -11 Torr. The measurements were carried out 
at T = 107 K which is below the pseudogap temperature 
T* . The pseudogap temperature of the sample may be 
determined from the temperature dependence of the re- 
sistivity shown in Fig. [5] From the deviation from the T 
linearity of the resistivity, it is approximately T* w 160 
K in agreement with the independent estimate from dif- 
ferent group 

By the Lorentzian fitting of the ARPES momentum 
distribution curve (MDC) shown in Fig. 21 we extract 
the self-energy £(0,w) (shown in Fig. [5] below) using a 
realistic tight-binding band dispersion (Eq. (0} below). 
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We confirm that the normal state self-energy exhibits an 
in-plane momentum anisotropy in both the elastic and 
inelastic parts. Details will be presented in Section II. 
From the obtained self-energies, we extract the function 
a 2 F by inverting the Eliashberg equation. We employ 
the adaptive maximum entropy method. Its formula- 
tion is given in Section III. We show that in spite of the 
anisotropy in E(#, oj) the a 2 F(9,u>) obtained from data 
from the cuts at different 0's collapse onto a single curve 
(Fig. [7]) at low energies with a cut-off at about 0.4 eV 
around the nodal direction changing smoothly to approx- 
imately 0.2 eV at 25 degrees. The variation of the cut-off 
is simply accounted for by the variation in the position 
of the bottom of the band with respect to the Fermi- 
energy with 9 with the intrinsic cut-off of the spectrum 
of about 0.4 eV. Some of these results were anticipated by 
an approximate calculation 2 ^ using an assumed spectra 
by analysis of ARPES in La2- a; Sr 2 ;Cu04, rather than the 
more reliable inversion method used here. These results 
put strong constraints on the momentum dependence of 
the fluctuation spectra, which we will discuss later. An 
implication of our results is that the anisotropy of the 
transport and spectroscopic properties originate simply 
from the anisotropy of the Fermi surface. The detailed 
results are presented in Section IV. In Section V, we will 
conclude with a discussion of the results, summary and 
remarks motivating further work. 



II. DEDUCING THE SELF-ENERGY 

The ARPES intensity, within the sudden approxima- 
tion, is given by 

J(k,w) = |M(k,i/)| 2 /M [A(k,w) + B(k,w)] , (2) 

where M(k, v) is the matrix element, v the energy of in- 
cident photon, f(u>) the Fermi distribution function, A 
the qp spectral function, and B is the background. We 
write the in-plane momentum k with the fcj_ perpendic- 
ular to the Fermi surface (FS) and the angle 9 measured 
from the nodal cut as indicated in Fig. [TJ Since E has a 
much weaker dependence on k± than u (which is verified 
in the experiments through the Lorentzian distribution 
of the spectra as a function of (k± — kp(9)), the spectral 
function may be rewritten as 



1 



£2(0,0;) 



7r[w-£(k)-Ei(0,w)] a + [£ 2 (0,w)] s 



A(0,k±,u) 

Here, £(k) is the bare dispersion, and Ei and E2 are, re- 
spectively, the real and imaginary parts of the self-energy. 
Along a cut with a fixed tilt angle 9, then, the depen- 
dence on k± is solely through the bare dispersion energy 
£(k). A has a Lorentzian form with respect to £(k). The 
self-energy can be directly extracted from the ARPES 
intensity by a Lorentzian fit if £(k) is known. 

We first determined the tight-binding (TB) dispersion 
of £(k) by matching the experimental Fermi surface with 



the four parameter dispersion. 

£(k x , k y ) = —2t(cosk x a + cosk y a) + At' cosk x acosk y a 

— 2t" (cos2k x a + cos2k y a) — //, (4) 

where a — 3.82 A is the lattice constant and fi is the 
chemical potential. We took t = 0.395, t' = 0.084, 
t" = 0.042, and /z = —0.43 eV, which are consistent 
with Kordyuk et al£2- Note that we neglected the bilayer 
splitting present in the Bi2212 compounds. At the pho- 
ton energy used (hp — 6.994 eV) in the laser ARPES, 
only the antibonding bands are observed and the bond- 
ing band is completely suppressed. The experimentally 
determined FS in comparison with that from Eq. Q is 
shown in Fig. [TJ The 6 cuts with the tilt angles 9 with 
respect to the (77, it) are also shown with the solid lines. 
To study the importance of the bare dispersion £ we also 
used the linear dispersion (LD) for comparison. 



£(k)=u,(0) [k x -k F (9)} 



(5) 



where Dp and kp were calculated from the tight binding 
dispersion of Eq. (|4]). 

In order are some technical points about the cuts per- 
pendicular to the Fermi surface. In the ARPES experi- 
ment setup, two angles are the control parameters: one 
is the tilt angle 9 which closely corresponds to the actual 
tilt angle 9 and the other is <f> which determine the k± 
such that 



ko = — ^ kl — , k x = kasmc/), k y — ko cos (j) sin 9 , (6) 

where Ekin is the kinetic energy of the photoelectons 
given by 



hv-W 



(7) 



where W = 4.3 eV is the work function. The k x and 
k y are the components of the wave-vector of the pho- 
toelectron with respect to the diagonal cuts. The k±, 
which is the distance from the (tt, tt) point can be simply 
calculated by the trigonometry. The <fi is varied by ap- 
proximately 30 degrees and the corresponding k± which 
go through the Fermi surface is shown by the thick bars 
in Fig. [T] A consequence of this experimental setup is 
that the actual tilt angle 9 is not constant as given by 
the 9. 9 = 9 for 9 = 0, but 9 begins to deviate from 9 as 
4> is varied when 9^0. The actual cuts therefore are not 
straight lines but slightly curved for large k± as shown in 
Fig. [l] but for k± near the Fermi surface the cuts point 
to the (tt, tt) direction. We disregard the difference be- 
tween the 9 and 9, and use 9 to denote the tilt angle in 
this paper. 

Typical ARPES intensity as a function of k± for fixed 
tilt angles and binding energy, referred to as the momen- 
tum distribution curve, and the dispersion as a function 
of the binding energy w are shown in Fig. |3J Each curve 
is shifted down for clarity as the binding energy is var- 
ied. The left and right plots are, respectively, for the tilt 
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FIG. 3: The representative momentum distribution curves (MDC) as a function of the momentum perpendicular to the FS, 
k±a/n, and their dispersion as the binding energy is varied. Each curve is shifted down for clarity. The left plot is for the tilt 
angle 5 and the right for 20 degrees. From the top to bottom are the MDC at the energy ui — 0.0005 eV with the step of 0.005 
eV up to 0.3455 eV. The qp coherence peak becomes suppressed as the energy is increased away from the Fermi level or the 
tilt angle increased from the nodal cut. 



angle 5 and 20 degrees. From the top to bottom are the 
MDC at the energy u> = 0.0005 eV with the step of 0.005 
eV up to 0.3455 eV. Note that the qp coherence peak 
becomes suppressed as the energy is increased away from 
the Fermi level or the tilt angle increased from the nodal 
cut. These MDC were fitted by equating the ARPES 
intensity with the spectral function given by Eq. ([3]) as 

I(6,k±,u)) = CA(9,k±,u)) + B(6,u) (8) 

to extract the self-energy as a function of the energy for 
a given cut. 

In Fig.H we show the MDC fits for the nodal (0 = 0), 
9 = 15, and = 25 degrees cuts and for the binding en- 
ergy u) = 0.0005 (the Fermi level) and uj = 0.0975 eV as 
representative cases. The shaded dots are the experimen- 
tal ARPES intensities and the solid dots are those used 
in the fitting. The red and blue lines are the spectral 
function of Eq. Q with the TB and LD bands, respec- 
tively. From the peak position and the width of the peak 
together with the bare dispersion £(k), we can determine 
the real and imaginary parts of the self-energy. 

The self-energies, — T,i(9,uj) and Y,2(9,ui) for u> > 0, 



determined this way at T = 107 K are shown in Fig. [5] 
The plot (a) was obtained using the TB and the plot (b) 
using the LD band. The real parts of the self-energy, 
Ei(w), cross the zero at progressively smaller energies 
as the tilt angle is increased. This feature is more pro- 
nounced in the TB dispersion analysis as can be seen by 
comparing the plots (a) and (b), and is better described 
by the TB because the band bottom can not be captured 
by the LD. The imaginary parts of the self-energy de- 
crease monotonically as w is increased up to u — 0.45 
eV. The elastic part, Y,2(9,u> = 0), clearly changes as 6 
is changed. The elastic qp scattering rate is momentum 
anisotropic in accord with previous works J^r— The func- 
tional dependence of the qp scattering rates on the 9 and 
lu for small uj may be analyzed in analogy with the Eq. 
([1]). We confirm that both elastic and inelastic qp scat- 
tering rates exhibit an anisotropy as a function of angle 
9 around the Fermi-surface. 

The extracted self-energies were used as an input to 
dudece the fluctuation spectral functions by inverting the 
Eliashberg equation. It will be presented in the following 
section. Before the detailed analysis, we note that the 
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FIG. 4: The representative momentum distribution curves as a function of the momentum perpendicular to the FS, k±a/n, 
at the tilt angles 6 — 0, 15, and 25 degrees, and the binding energy ui — 0.0005 and 0.0975 eV. The shaded circles are the 
experimental ARPES intensities and the solids ones are those used in the fitting. The red and blue lines are the fitting with the 
TB and LD bands, respectively. Note that the MDC is not Lorentzian when plotted as a function of k± with the TB dispersion 
as can clearly be seen from (f). 



relation 



dT, 2 {e,uj) 
duo 



ira 2 F(9,uj) 



(9) 



holds for uj ^> T, as may be deduced from Eq. (112]) be- 
low. The extracted £2(6', a;) from the TB collapse onto a 
single curve up to uj < 0.2 eV as can be seen from Fig. [5] 
This, together with Eq. ([9]), suggests that the Eliashberg 
function a 2 F(9, ui) would yield a single curve indepen- 
dent of the tilt angle 9 for T < ui < 0.2 eV. The detailed 
analysis presented in the following sections will establish 
that the Eliashberg function a 2 F(9, uj) at different angle 
9 collapse onto a single function of uj with the intrinsic 
cut-off of about 0.4 eV or the energy of the band bottom 
in direction 9 with respect to the Fermi energy, whichever 
is smaller. This is our key result which we will turn to. 



III. DEDUCING THE FLUCTUATION 
SPECTRA: FORMULATION 

Neglecting vertex corrections (justification for this 
will be given in the concluding section), the self-energy 
S(k, uj) may be written as 



/OO POO 
de / de 
-OO J — OO 



, f(e) + n(-e') 



e + e' — uj — 18 
^A(k',e)a 2 (k,k')F(k',k,e'), (10) 



where / and n represent the Fermi and Bose distribu- 
tion function, respectively. a(k, k') is the matrix element 
for scattering fermions with the fluctuations of spectral 
weight F(k',k, e). After the integral over k' ± , the real 
and imaginary parts are given by 



Ei(0,w) 



,// / de m+n{ - e 'K*F(e,e>),(ll) 



e + e — ui 



Z 2 (9,uj) = 7t / de' [f{uj - e') + n{~e')]a 2 F(9,e')ll2) 



where the Eliashberg function, or the bosonic coupling 
spectrum, a 2 F(9, e') is given by 



a 2 F(9,e') 



The Vf{9') is the angle dependent Fermi velocity and the 
bracket implies the angular average. We can extract the 
"averaged" spectrum a 2 F(9, e') by inverting either Eqs. 
(fTTl) or (|T2"j) . Thus information only about the angle de- 
pendence of the average over the product of the squared 
matrix element and spectra as constrained by energy- 
momentum conservation are obtained. Both equations 
must give the same results for a 2 F(9,ui) provided the 
real and imaginary parts of the self-energy satisfy the 
Kramers-Kronig (KK) relation. In the present case, how- 
ever, the extracted spectra do depend on whether it is ex- 
tracted from the real or imaginary part. It implies that 
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FIG. 5: The real and imaginary part self-energies at T — 107 
K for the tilt angles 9 = 0, 5, 10, 15, 20, and 25 degrees. -Si 
and E2 as a function of the positive energy ui are shown. Plot 
(a) shows the extracted self-energy using the TB dispersion, 
and (b) shows, for comparison, one using the linear dispersion. 



the extracted real and imaginary parts of the self-energy 
do not satisfy the KK to the required accuracy. It turned 
out that the a 2 F from the real part self-energy is more 
reliable as shown below, probably because the peak posi- 
tion of MDC is better determined than the width. This 
point was also noted by Shi et aL- 

A few techniques have been devised to invert the 
Eliashberg equation, (fTTT) or (fT2|) . We employ the max- 
imum entropy method (MEM) i 18 ' 29 It is a useful tech- 
nique to overcome the numerical instability in the direct 
inversion, by incorporating the physical constraints into 
the fitting process. The MEM minimizes the functional: 



»=i 
p 00 

S = de' 
Jo 



07 



a 2 F(e') - m(e') - a z F{e')\n 



m(e') 



(15) 



(16) 



where Njj is the total number of data points, is the 
experimental data for the real part of the self-energy at 
the energy Wj, w») is defined by Eq. (fTTj) . and Uj 

is the error bars of the data. The entropy term imposes 
physical constraints to the fitting and is maximized when 
a 2 F(e') — rale 1 ), where m(e') is the constraint function. 
To minimize Eq. (fPlj) for a given m(e'), the Eliashberg 
function a 2 F(uij) is optimized in each iterative step by 
updating it as 



5F ]= -Y,A 



6L 

3k W k 



A jk = 



5 2 L 
6F j 5F k 



(17) 



where Fj — a F(uij). The matrix inversion of A~ k was 
performed using the singular value decomposition tech- 
nique. Details of the algorithms can be found in Ref^S. 

In usual MEM approaches, the constraint function is 
chosen based on our best a priori knowledge for the spe- 
cific system and remains unaltered. The obtained results 
for a 2 F(e') do depend on how the m was chosen. In 
order to decrease the dependence on the constraint func- 
tion and also to better represent the physical system, we 
employ the adaptive MEM which was implemented as fol- 
lows: For a chosen constraint function m(e'), the iterative 
minimization is performed as described above. After the 
convergence, the m(e') is updated to a linear combination 
of m(e') and a 2 F(e') to reflect the nature of the solution. 
For this new constraint function, the minimization is per- 
formed again via iterations. This second optimization is 
repeated until the convergence is reached. The double it- 
erative adaptive MEM decreases the dependence on the 
choice of the constraint function and improves the over- 
all fitting quality. The quality of the adaptive MEM can 
be seen by comparing the real part self-energy from the 
ARPES data and the self-energy expression of Eq. (fTT|) 
as shown in Figs. [6l 



The x 2 is the error and S is the generalized Shannon- 
Jaynes entropy defined below. The multiplier a is a de- 
terminative parameter that controls how close the fit- 
ting should follow the data while not violating the phys- 
ical constraints. When a is small, the fitting will follow 
the data as closely as possible at the expense of a noisy 
and/or negative Eliashberg function, and when a is large, 
the extracted Eliashberg function will not deviate much 
from the constraint function m(e'). For a given tilt angle 



IV. FLUCTUATIONS SPECTRA 

As noted previously, the real part of the self-energy is 
determined more accurately from the ARPES intensity, 
and as the tilt angle is varied the tight-binding dispersion 
is more reliable. We therefore present in Fig.[7]the Eliash- 
berg function a 2 F(9, ui) extracted from the real part self- 
energy determined using the TB dispersion shown in Fig. 
[5] Note that structures in £(6*,oj) are reflected in the 
coupled fluctuation spectra a 2 F(0,ui). The Eliashberg 
functions increase approximately linearly from zero as u> 
is increased and have a peak at w ~ 0.05 eV, flattens 
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FIG. 8: The angle dependence of the Fermi velocity and band 
bottom calculated from Eq. @. 



FIG. 6: Comparison between the real part of the self-energy 
and the MEM fitting of Eq. (TT} for the tilt angles = and 
15 degrees. The squares are the extracted real part self-energy 
of Fig. Oa), and the solid lines are the MEM fitting. 
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FIG. 7: The Eliashberg function extracted from the real 
part of the self-energy at T = 107 K. Notice the remarkable 
collapse of the Eliashberg functions below uj ~ 0.2 eV for 
different tilt angles. 



peak energy remains unaltered, but the cut-off energy 
decreases monotonically. At 8 — 25 degrees, uj c « 0.25 
eV. With Laser ARPES the data at large angles is not 
obtainable due to kinetic constraints; moreover the prob- 
lems presented in the analysis due to interference of data 
due to the contributions from the next zone, as is ob- 
vious from Fig. [TJ become more serious. An important 
point may be noted from comparison of Fig. [7] and Fig. 
[8] where the position of the band bottom with respect 
to the Fermi-energy is given. The latter is about 1 eV 
in the nodal direction and about 0.25 eV at 9 25 de- 
grees. As was shown in a simple calculation^ and as 
is natural, the band-bottom serves as an effective cut- 
off in the fluctuation spectrum when it is lower than the 
intrinsic cut-off. So one important results is that the in- 
trinsic cut-off of the spectrum is about 0.4 eV. An even 
more important result contained in Fig. [7] is that quite re- 
markably the fluctuation spectra for different tilt angles 
collapse onto a single curve below the angle dependent ef- 
fective cut-off. Note that the corresponding self-energies 
do change even at low energies as the angle is varied as 
can be seen from Fig. From Eq. (jTTj) it follows that a 
nearly isotropic a 2 F{9, uj) produces an angle dependent 
self-energy purely from the effects of the dispersion of the 
bare band. 



above 0.1 eV, and vanishes above a cut-off energy uj c 
which does depend on 9. The noisy feature above 0.35 
eV is the numerical artifacts of the singular value decom- 
position and MEM. The cut-off energy is approximately 
0.37 eV along the nodal cut. These are in good agreement 
with Schachinger and Carbottei^ who analyzed the nodal 
cut Bi2212 laser ARPES data for the same sampled 

Now we discuss the new results of the analysis of the 
off-nodal cuts. The Eliashberg functions obtained are 
presented in Fig. [7] for the tilt angles 9 — 0, 5, 10, 15, 
20, and 25 degrees. As the tilt angle is increased, the 



The Eliashberg functions may also be extracted from 
the imaginary parts of the self-energy using Eq. (|12p . 
The results are given in Fig. |H] for the tilt angles and 
15 degrees for representative cases. As alluded before, 
they do not agree exactly with those from the real parts 
of the self-energy. Especially the features around the 
cut-off are exaggerated in the a 2 F(8,ui) from the imagi- 
nary part self-energies. The a 2 F(8,uj) do not collapse as 

neatly as those from the real parts. Also shown are the 
_i OSaM f th tnt i Th _i8M£) d 

7T OUJ ° 7T OUJ 

a 2 F{uj) determined from £2(0;) agree overall as expected. 
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FIG. 9: The Eliashberg function deduced from the imaginary 
parts of the self-energy. Results for 8 = and 15 degrees 
are given for representative cases. Also shown are 
They agree overall as expected. 
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V. DISCUSSION OF RESULTS, CONCLUSIONS 
AND GUIDES TO FURTHER EXPERIMENTS 

In this section, we discuss the implications of our find- 
ings which point to future experiments and analysis. 
They are the angle independence of the Eliashberg func- 
tion a 2 F(6,uj), i.e the angle averaged product of the ma- 
trix element and the spectral function of the fluctuations, 
as defined in Eq. (fT3|) , the possible physics of the low en- 
ergy bump around s» 0.05 eV, the vertex corrections to 
the Eliashberg equation, and various other assorted is- 
sues. The finding that the Eliashberg function is angle 
independent below the cut-off uj c puts an important con- 
straint on the microscopic understanding of the cuprates. 
For example, it can put a limit on the correlation length 
£ for the commonly assumed form of the antiferromag- 
netic (AF) fluctuations. A phenomenological form of the 
overdamped AF fluctuations may be written as 2 -^ 



XAF(k,k',w) = 



a£ 2 lv/ujai 



((k-k'-Q) 2 e 2 + l) 2 + (^/ W AF) 



r(18) 



Xaf{9, <j->) can be obtained after the integral over 9' with 
both k and k' on the Fermi surface, that is, 



(19) 



where k and k' have the angle 9 and 9' with respect to the 
nodal cut, respectively. A straightforward calculation re- 
veals that a weak 9 dependence of xaf{0,u) means that 
£/a <C 1 for uj < ujaf, where ujaf is the characteristic 
AF energy scale. This is shown in Fig. [TU]for £/a = 1/tt 
with the blue and for £/a = 1 with the red lines. For each 
£, the angles are = 0,5, 10, 15, 20, and 25 degrees from 




FIG. 10: The model Eliashberg function calculated from the 
overdamped AF fluctuations of Eq. (|18[) . The red and blue 
lines are for the AF correlation length {/a = 1 and re- 
spectively. For each f, the angles are from above 9 = 0, 5, 
10, 15, 20, and 25 degrees. As £ is decreased the fluctuation 
spectra become momentum isotropic. 



above. As expected, xaf{0,uS) becomes angle indepen- 
dent as £ is decreased. The collapse of a 2 F(9,ui) implies 
that if it is due to the AF fluctuations, the AF correla- 
tion length should be £/a < 1/tt for Bi2212 at T — 107. 
This is an upper limit since we have not included the an- 
gle dependence of the matrix elements of coupling of AF 
fluctuation to fermions which enters Eq. (fl~3|) . A correla- 
tion length less than a lattice constant is indicative either 
that there is negligible spectral weight in AF fluctuations 
or that they couple to single-particle excitations weakly 
enough to have no observable effect. This implies among 
others that the antiferromagnetic fluctuations may not 
be underlying physics of the deduced fluctuation spec- 
trum. This also means that the scenario of "hot" and 
"cold"- spots on the Fermi-surface is inapplicable to the 
cuprates. 

An angle-independent o 2 F(9,oj) is consistent with the 
quantum-critical spectra hypothesized in the marginal 
Fermi-liquid description of cuprates and recently derived 
microscopically 31,32 to be the spectra in the quantum- 
critical region of the phase diagram of the cuprates 
due to the quantum melting of the loop-current order 
observed 3 ^— in the underdoped region of the cuprates. 
There is one aspect of the deduced a 2 F(9,uj) which is 
not given by the theory. This is the low energy bump at 
the energy of u) « 0.05 eV. The presence of this bump 
may be seen directly in the deduced Im E(0, u)) which is 
not linear as a function of ui for low uj 3> T. The linearity 
may be expected on the basis of such theory and earlier 
ARPES experiments which, however, do not have the 
high resolution of the present experiments. If the bump 
occurs only in samples studied in the pseudogap region, 
there must exist collective modes special to this region of 
the phase diagram. This can be checked by equally high 
resolution data in samples in the quantum-critical region 
or the overdoped region of the phase diagram. 

We now discuss the validity of the neglect of vertex cor- 
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rection in the relation of the self-energy and the fluctu- 
ations. The exact relation between the vertex correction 
<5A(k, q, uj, v) and the self-energy is that 

SA(k,q,w,z/> = E(k + q,w + i/) - E(k,w). (20) 

As discussed, the experiments give at energies up to 
about 0.2 eV, that the self-energy has a momentum de- 
pendence due only to the anisotropy of the Fermi veloc- 
ity, which from Fig. [5] is only about 20 %. This may 
be regarded as negligible. It therefore follows that the 
momentum dependence of the vertex correction is negli- 
gible. This means that all the essential conclusions about 
the momentum independence of the fluctuation spectra 
arrived at here remain valid. As regards the frequency 
dependence, there is indeed a vertex correction which in 
dimensionless form is of 0(a 2 uj c /W), where uj c is the 
upper-cutoff of the fluctuations and W is the electronic 
bandwidth of about 2 eV. For a 2 ~ 1, there are then 
vertex corrections of 0(1/5). This means that our quan- 
titative conclusions have a validity no better than about 
20 %. 

One final remark pertains to the information that the 
deduced spectra offers for the fluctuations spectra re- 
sponsible for the superconductive instability. The nor- 
mal state single-particle self-energy has the full symme- 
try of the lattice. It then follows, as is visible from Eq. 
(jTQJ) that the u £ = 2" part of the fluctuation spectrum, 
which determines the symmetry of the superconductiv- 
ity in cuprates is not visible through study of the nor- 
mal state. It is unlikely that a completely different form 
of the spectra is responsible for the normal state and 
the superconducting state. A test of this conjecture and 
the deduction unambiguously of the spectra responsible 
for superconductivity requires a study of the variation 
of the single-particle self-energy as a function of energy 
on going from the normal to the superconducting state. 
This would be a generalization^ of the McMillan-Rowell 



procedure^ to d-wave superconductors requiring high 
resolution ARPES. Such work is in progress. 

To summarize, we have presented the self-energy anal- 
ysis from the ultra high resolution laser ARPES on the 
slightly underdoped Bi2212 high temperature supercon- 
ductors. The self-energy was determined along the nodal 
and off-nodal cuts in the normal state. Both the elastic 
and inelastic quasi-particle scattering rates exhibit in- 
plane momentum anisotropy. The deduced self-energies 
were then used as an experimental input to invert the 
Eliashberg equation to extract the product of the fluc- 
tuation spectra and the coupling to the single parti- 
cles. The adaptive MEM was used. The high resolution 
ARPES data together with the realistic tight-binding 
band dispersion enabled us to determine high quality self- 
energy and the Eliashberg function. At the temperature 
T = 107 K between the superconducting T c and pseudo- 
gap temperature T* , the Eliashberg functions for differ- 
ent tilt angles collapse onto a single curve up to the up- 
per cut-off despite the angle dependent self-energy. The 
cut-off has an intrinsic value of about 0.4 eV or the band- 
bottom energy with respect to the Fermi level if it is less 
than about 0.4 eV. The implications of our results have 
also been discussed. 
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